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Abstract 

Assuming a simple spherical relativistic mean field model of the nucleus, 
we estimate the width of the antiproton-neutron annihilation (F^) and the 
width of antiproton-proton (Tp) annihilation, in an antiprotonic atom system. 
This allows us to determine the halo factor /, which is then discussed in 
the context of experimental data obtained in measurements recently done on 
LEAR utility at CERN. Another quantity which characterizes the deviation 
of the average nuclear densities ratio from the corresponding ratio of the 
homogeneous densities is introduced too. It was shown that it is also a good 
indicator of the neutron halo. The results are compared to experimental 
data as well as to the data of the simple liquid droplet model of the nuclear 
densities. The single particle structure of the nuclear density tail is discusssed 
also. 
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I. INTRODUCTION 



The LEAR measurements at CERN |jT|,0], show large neutron halo factors for selected 
nuclei. The observations of this factor were performed by counting those nuclei which were 
born in annihilation reactions in antiprotonic atoms. The following processes were observed: 

p +i Xn — ^^z^ Xn-1 + {it} . 

Here, zXn denotes the nucleus with the atomic number Z, the mass number A and the 
number of neutrons A^. The annihilation products consist of a few (usually 4-5) pions {vr}. 

In reactions shown the resulting nuclei have the number of neutrons or the number of 
protons diminished by one as compared to the mother nucleus. Counting the annihilation 
products resulting from nuclear neutrons Afn and from nuclear protons Afp, one can calculate 
the following halo factor / = {Z / N)N'n/ Mp. For the case of "cold" annihilation reactions 
which take place at the peripheral region of the nucleus and do not excite the final nuclei, 
the factor / is supposed to depend strongly on the ratio of neutron and proton densities. 

In order to describe the phenomenon of annihilation of the antiprotonic atomic state on 
the mother nucleus we have done the calculation of the halo factor / and the density ratios 
deviation which we shall define later. As a basis of our calculation serves the relativistic 
mean field model in which the field equations are solved numerically. These solutions are 
exact even far outside the nuclear surface. This precision is critical for the calculation of the 
quantities of interest. Both quantities / and have been considered for nuclei: ^^Ni, ^^Zr, 
96 Ru, i30Te, i^^Sm, ^^^Sm, ^^^Yb, ^ssjh, 238u. All these nuclei were studied in Ref. [§0. The 
halo factor / shows also some dependence on the separation energies of the last nucleons 
(proton or neutron). 

This paper is organized as follows. In the next section we define and discuss the model 
independent functions / and h which are supposed to characterize halo phenomena. An 
example, the droplet model calculation of h is shown. The third section contains a short 
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presentation of the relativistic mean field model which we have used to determine / and h 
factors. 

In the Summary of the paper we show and comment on the results of our calculations 



II. BASIC DEFINITIONS 



In this section we define basic quantities which are suitable in discussing halo phenomena. 
The first subsection collects the approximations made in halo factor definitions. The second 
subsection introduces new quantities. These are average densities calculated from the density 
distribution and the deviation of average densities ratio from the homogeneous neutron and 
proton densities ratio. An example of this is shown and reviewed. 



The simple expression for the antiproton absorption width for the antiproton which 
occupies a state (s) of an atom is 



The subscripts n and p are for neutrons and protons respectively. The wave function ^'^(r) 
of antiproton is determined in the Schrodinger model of hydrogen-like antiprotonic atom 
within a point nucleus. The approximation of the point like nucleus is motivated by the 
fact of a large mean distance B of the antiproton from the nucleus. Taking the mass ratio 
of the antiproton to the electron mp/me ~ 2000, the main (supposed) quantum number 
n of antiproton in the atom eg., n = 8, and the charge of the nucleus in question, is 
B K, i?o/2000Z, where i?o, is the Bohr radius of the electron in the hydrogen atom. For 
the cases which we are considering here this gives an average antiproton-nucleus distance of 
the value S 30 — 40fm. As it is seen, this distance exceeds a few times past the nuclear 
radius R ~ A^/^fm. This motivates our point nucleus approximation. It is possible to 
solve the Dirac or the Schrodinger equation for the more complicated case of finite or even 



A. Halo factor 




(1) 
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deformed nucleus. However, this will not change significantly the relative measures which 
we shall introduce and apply in this present paper. 

The factor P(r), in the definition of / describes pion escape probability plus other effects 
and it will be discussed later. Here, we may say only, that in the present model it does not 
infiuence the calculation and with very high accuracy it may be assumed as P{r) ^ 1. 

We define the halo factor as |^ 

r Z EsK 

A similar definition was used in [|l],|2|. The summation run through all the antiprotonic 
atomic states for which one expects the large annihilation widths. Usually it is assumed 
that only one state with {n, I) ranging from (6, 5) for Ni nucleus to (9, 8) to U nucleus 
contributes to these widths [§| . In the present calculation we sum a few (three to five) states 
only. The summation stops if the addition of the successive term does not change the 
value of the sum. 

It is worthwhile to mention the weak dependence of our / definition on the structure 
of the antiproton-nucleon (pN) interaction. It is assumed that both proton and neutron 

widths depend on the imaginary part of the complex optical potential W(r) which is 
responsible for the absorption of the antiproton. We use the approximation that the optical 
potential is proportional to the density of nucleons: W{r) = ap{r). Here the parameter a 
is the pN interaction length and is taken as a = + ia^, where both real and imaginary 
parts of a are constants ||^. In this way we obtain 

Tn{p) ~ an{p) p{r) . (3) 

Because both and in the definition of halo factor / (Eq. |^) appear in a ratio and 
thus depend weakly on the specific antiprotonic atomic state one concludes that the optical 
potential dependence of / is washed out. In this sense / is the pN interaction independent 
measure of the annihilation ratio and describes rather the geometrical properties of the 
nucleon density distributions. 
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In fact, the definition of / consists of an/^m ap which measures the cross sections an- 
nihilation ratio a{pn) / (j[pp) for the antiproton annihilation on nucleons in a nuclear medium. 
The value of this ratio for a given nucleus is very hard to estimate in an experiment or on 
theoretical grounds. In an early paper on nucleon-antinucleon optical potential by Bryan 
and Philips ||^, the relative antiproton p-capture rates on neutrons and protons in singlet 
and triplet states was estimated from measurements of p annihilation at rest in hydrogen 
and deuterium. From these one can calculate the ratio of capture on neutrons to the capture 
on protons. In singlet states it is close to 0.6 whereas in triplet states it exceeds 1. The 
measurements by Bugg et al. show that this ratio is close to 0.65 in case of the nucleus 
^^C. Other authors give us an imaginary potential ratio the value ~ 80%. In the paper 
by Kalogeropoulos and Tzanakos 0, one can find a similar value a{pn)/a{pp) = 3/4 as 
measured in deuterium. These different numbers suggest that the annihilation ratio may 
change from one nucleus to the other and in our opinion, it is not possible to scale the results 
using these numbers - the ratio may show a dependence on proton and neutron numbers Z 
and N as well as the state in which the annihilation takes place. 

Therefore, in the calculations which follow we do not scale our results. As we shall see 
the results show some similarity to the experimental data. This supports our choice of the 
RMF model calculation and shows that the ratio Qmanf^Tnap is close to unity for our 
model. 

To select the peripheral annihilation acts (which one observes in pX experiment) one 
should correct the definition of the halo factor on the pion escaping probability Pn,esc{f) and 
the deep hole creation probability Pfih{r). One has to take into account only this part of 
the annihilation width F which is responsible for the cold annihilation — the annihilation in 
which pions escape to infinity without exciting the rest of the system. On pure geometrical 
grounds this is shown to be a power function of r which agrees roughly with very accurate 
calculations presented in Ref. [^]. Assuming the dependence Pn,esc ~ r'', where k > 0, 
one can see from the Eq. |^ that the summation of widths F will run over nearly the same 
spectrum of powers of r as without this factor. It is therefore allowed to stay with the old 
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formulas in which P{r) is constant. 

Another factor which probably enters the expression for F is a deep hole creation factor 
Pdhi'f') 0- Its role is again the preservation of cold acts of annihilation. The physics behind 
this is the following. A number of annihilation acts may occur on a deeply situated nucleon 
levels. This leads to a rearrangment of the nucleonic orbits and it may happen that the final 
system will show very small or positive(!) binding energies of the last nucleons. After its 
emission, the nucleus in question goes out of the observation range and can not be counted 
as the peripheral halo product. The dependence of the function P^h on r seems to satisfy 
the power law analogous to pion escaping probability Pn,esc ~ r'', and as it was said before, 
we do not introduce this into our calculation. These are all the approximations which we 
have assumed in our paper while evaluating the halo factor. 



B. Average densities 

Consider an average density p defined by the integral 

P=^ I dTp\r), (4) 

where p(r) is the sphericaly symmetric nucleon density distribution and dr in the integration 
volume element. N is the number of particles in the system and is given by 



N = JdTp{r), (5) 
In the case of pure Fermi distribution 

p(r) = po 1 + exp , (6) 



the quantity p can be evaluated approximately with assumed accuracy. In the above formula, 
a is the width of the nuclear surface, R is the radius and po is the central density of the 
nucleus. The expression valid to second order in expansion parameter a/R can be obtained 
from the Elton's formula (pages 106-107 of Ref. []TU|)- To use this formula one has to modify 



the expression for the average density to the following form 



pr^ F2{k)-2F^{k), (7) 



where k = R/a and Fn{k) the Elton's integral is defined by 



JO (1 + exp(x — k)) 



From the formula mentioned before, one has to second order m. a/R 



p^Po(l-3- + y(-n. (9) 

A relatively large, linear term —3a/R, which appears in this expression makes the p a sen- 
sitive indicator of the peripheral properties of the nuclear densities. 

On the basis of both average densities p„ and Pp one can define the deviation of their 
ratio from the homogeneous densities ratio. This is 

h = l-^^^. (10) 

N Pp ^ ^ 

Inserting into the last equation shows 

/i~3^!i^, (11) 

where a„ and Op are the diffusness parameters of the neutron and proton distributions 
respectively and we have assumed Rn ^ Rp = R- 

Figure |I] shows an example of the deviation h as calculated for all of the viewed nuclei 
in the case of the liquid droplet model densities |Tl|. The deviation h was calculated from 
the parameters of this model and, the geometric dependencies of the a and R parameters 
entering the Fermi density distribution (see Eq. P) and the corresponding parameters of 
the liquid drop distribution [jl2|-Q. Figure ^ shows two kinds of data. The impulses show 



exact (full lines) and approximate (dashed lines) results respectively. One can see the good 
agreement of both exact, and approximate, calculated from Eq. (p!0| ) of the density deviation 
h values. Large positive values of h can be seen in the case of neutron rich nuclei like Te, 
Yb, Th and U. The deviation h is small and positive in cases of neutron deficient nuclei: 
Ni, Zr, Ru. Latter we shall see a similar behaviour of h for more realistic calculations done 
with the RMF model (see Figure 



III. RELATIVISTIC MEAN FIELD MODEL 



In this section we describe very shortly the relativistic mean field (RMF) model which we 
have used in the calculations of the neutron and proton densities entering both halo factor 
/ and the average density deviation. 

The RMF theory starts from a lagrangian consisting of nucleonic and mesonic degrees 
of freedom In some sense it seems to be more fundamental than other mean field 



models like Skyrme-Hartree-Fock and Gogny-Hartree-Fock models. It gives the relativistic 
treatment of nucleonic and mesonic variables and a proper description of the spin-orbit 
interactions. Nevertheless it is still an effective phenomenological method based on the local 
densities and fields. 

The RMF theory was successfully used to reproduce parameters of the nuclear matter 
and some properties of finite nuclei like binding energies, mean square charge radii and 
quadrupole moments |lT6|jr7| . 

This theory is based on the following field lagrangian density |[T5|JT6|JT8| -pO 



—Q^'Tl^^ + ^ml uj^^Uf, - ^ oj^, (12) 
-^R'"'R^^ + p^p^ - Qp^i Yr^i Pf, 

The fields belong to nucleons (Dirac spinor field the low mass isovector- vector me- 
son {p^;Rp_u), isoscalar-vector {u^;Q^i^), scalar a and to the massless photon vector field 
{A^; F^y). The field tensors are given by 

= d^'u" - d^'u" , (13) 
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F"" = d'^A" - d^'A" . (15) 

The Dirac spinors ipi of the nucleon and the fields of a, p and u mesons are solutions of 
the coupled Dirac and Klein-Gordon equations which are obtained from Eq. |l^ by means 
of the classical variational principle and they are then solved by iteration for the case of 
the spherically-symmetric systems of nucleons. This iteration goes through the following 
steps: we start from an estimate of meson and electromagnetic fields and we solve the Dirac 
equation getting the spinors ipi- The spinors are then used to obtain the densities. The 
latter serve as source to solve the Klein Gordon equations and achieve the new estimations 
of the meson and electromagnetic fields. The parameters used in our calculation are the 



same as in Refs. [18,19,21 



IV. RESULTS 



On the basis of the relativistic mean field model which was presented in the previous 



section and is described in detail in |[r5| , p!8|J19| , |2l|j2Q , p!6|l we have calculated the nucleon den- 



sities entering the Eq. ^ for the absorption width and the halo factor Eq. |^. We show that 
some of the nuclei are good candidates for the neutron halo systems and have promising 
large logarithm of neutron to proton density ratio l{r). This ratio which we define as 

l{r) = log(p„(r)/pp(r)) . (16) 

is a sensitive quantity and it suits to relate densities especially in a peripheral nuclear region 
where both neutron and proton densities are very small. Figure shows this ratio for all 
nuclei under consideration. 

In Figure |^ one can see nuclei for which the logarithm /(r), for r > 2R, has the value 
of few orders. This means that the neutron density is there still greater then the proton 
density. This indicates a possible large halo factor /. 

One has to remember that the microscopic density p is the sum of the single particles 
densities over all occupied states 
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p(r) = Y: Pu{r) . (17) 

It is interesting to study the dependence on r of the contribution of the single particle 
orbitals v = {nlj) to the total density? Or in another words what is the collectivity of the 
nuclear density tail? The results of the RMF calculation for some nuclei in which the halo 
factor is discussed are presented in Figure ^ It is seen that for large distances of 10-14fm 
only a few valence orbitals contribute to the density tail. Even more, for some nuclei in 
which / is large, one neutron orbital exhausts 90% of the nuclear density. Only in the case 
of ^^^Sm in addition to the neutron orbital I/19/2 there appears 2^5/2 proton orbital with the 
comparable amplitude. In this nucleus one really observes very small neutron halo factor / 
(proton halo). 

We have to add from our numerical experience that the asymptotic behaviour of the 
single particle density depends dramatically on the single particle energy. For the orbitals 
less bounded the large r tail of density is longer. So the halo effect will not be the only 
probe of the surface width of the neutron and proton distributions but also a crucial test of 
the single particle structure foreseen by theoretical models. This also means the inclusion 
of a nuclear shape deformation into our model could be important. 

Figure § shows the dependence of the logarithm of the total density of baryons against 
the distance r measured from the center of the nucleus. The logarithm of the density in the 
peripheral area of the nucleus (r > 8) takes very small values and is nearly linear function 
of r with a steep negative slope. All the slopes for different nuclei are similar. 

We now consider another two quantities which characterize the real density distribution. 
The first is the deviation of density ratio h from the homogeneous density ratio ([T0|). In 
contrary to the case shown already in Figure |l|, this is calculated for real density distributions 
determined from the RMF model. In analogy to Eq. |10|, it is defined through the average 
densities calculated in RMF model. The second quantity we considered is the celebrated 
halo factor / (Eq. 0). Both quantities give the measure of the density distribution in the 
nucleus. These are shown in Figure |^ and Figure ^ respectively. 
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One can easily identify a halo nuclei. In Figure |^ you can see an interesting case of Ru 
nucleus for which the relative density deviation h is negative. This suggests the lack of the 
neutron halo for this nucleus. At the same time the halo factor / ~ 2. We concluded that 
the deviation h is rather a rough indication of the neutron halo in nuclei. 

The calculated halo factors / are shown in Figure ^ and in Table 1 where in addition we 
have also displayed the experimental values /exp- These values were extracted from Ref. 0. 
Except for extreme cases of ^^Ru, ^"^^Sm and ^''^Yb one can observe a good agreement 
between calculated factor /rmf and the experimental data. The case of ^"^^Yb it is a special 
one. It shows a too low theoretical halo factor 3.1 and a very high measured value of 8.1. 
This is the case discussed also in [Q. The case of samarium nucleus, ^^"^Sm, shows / which 
probably indicates the proton halo instead of neutron one. In our calculations we obtained 
in this case / = 1.5. The similar situation of large theoretical / is observed in the nucleus 
96Ru. 

V. SUMMARY 

In the relativistic mean field model we have calculated the halo factor / and the factor 
h - the average density ratio deviation from the ratio of homogeneous density distributions, 
for nuclei in which we observed the neutron halo. 

Both quantities / and h indicate the possible halo effects in most of the studied nuclei 
but in some cases (see eg., Yb nucleus) they differ from experimental data significantly. 
This fact shows that one has to extend the RMF nuclear density model. The extension 
may include eg., the deformation of the nuclear systems or /and the pairing interaction. The 
finite deformation in all considered cases may change the picture a little bit of the neutron 
halo. 

The new parameter characterizing the peripheral properties of the nuclear density dis- 
tributions, which we have called the average density ratio deviation /i, is proportional to the 
difference of the proton and the neutron surface diffuseness parameters. Like /, it points 
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out these nuclei which show the neutron halo properties. 

We have shown that the nuclear density tail manifests a single particle nature. Therefore, 
the neutron halo is the single particle effect and not a collective one. 

Some of our predictions, as compared to the experimental data, fail but most of them 
show correlations which are very promising ones. It is hard to explain the source of existing 
divergences. A possible explanation may be the lack of deformation in our model calcula- 
tions. It is also possible that the inclusion of pairing interaction to this theory may improve 
the results. Such calculations are in progress. 
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TABLES 



TABLE I. Experimental and calculated halo factors. 



Nucleus 


/cxp 


/rmf 


58Ni 


0.9 


1.2 


96Zr 


2.6 


2.3 


96Ru 


0.8 


2.3 


i30Te 


4.1 


3.5 




< 0.4 


1.5 




2.0 


3.0 


176Yb 




3.6 


232 rji]^ 


5.1 


5.5 


238 u 


6.0 


5.0 
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FIGURES 

FIG. 1. The deviation h of the neutron to proton density ratio for Fermi density distributions 
in the case of hquid droplet model. The impulses show exact (full line) and approximate (dashed 
line) results. The data corresponds to the following nuclei: ^^Ni, ^^Zr, ^^Ru, "'^^''Te, ""^^^Sm, ""^^^Sm, 

176Yb, 232Th, 238U. 

FIG. 2. The logarithm of the neutron to proton density ratio calculated on the basis of 
relativistic mean field model. The Z(r) curves are shifted on a value 6 given on the right hand side 
of each curve. 

FIG. 3. The partial density ratio for single particle neutron (n; full line) and proton (p; 
dashed line) orbitals. Arrows show the locis of the nuclear root mean square radii of the charge 
distributions. 

FIG. 4. The logarithm of the baryon density against the distance r (in fm) measured from 
the center of the nucleus. 

FIG. 5. The deviation h of the neutron to proton density ratio. The nuclei are the same as in 
Figure 2. 

FIG. 6. Halo factor /. The nuclei are the same as in Figure 2. The experimental data 
(squares) were taken from [^. 
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